Two-body correlation functions 
in nuclear matter with np condensate 
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The density, spin and isospin correlation functions in nuclear matter with a neutron-proton (np) 
condensate are calculated to study the possible signatures of the BEC-BCS crossover in the low- 
density region. It is shown that the criterion of the crossover (Phys. Rev. Lett. 95, 090402 (2005)), 
consisting in the change of the sign of the density correlation function at low momentum transfer, 
fails to describe correctly the density-driven BEC-BCS transition at finite isospin asymmetry or 
finite temperature. As an unambiguous signature of the BEC-BCS transition, there can be used the 
presence (BCS regime) or absence (BEC regime) of the singularity in the momentum distribution 
of the quasiparticle density of states. 
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Introduction. Neutron-proton pairing correlations 
play an important role in a number of contexts [l|, |3|, 
including the study of medium mass N ^i Z nuclei 
produced at the radioactive nuclear beam facilities H 
and the process of deuteron formation in medium-energy 
heavy ion collisions ^. For not too low densities, np 
pairing correlations crucially depend on the overlap be- 
tween neutron and proton Fermi surfaces and even small 
isospin asymmetry effectively destroys a condensate with 
np Cooper pairs due to the Pauli blocking effect [a, S 13 ■ 
However, under decreasing density, when neutrons and 
protons start to bind in deuterons and the spatial sep- 
aration between deuterons, and between deuterons and 
extra neutrons is large, the Pauli blocking loses its effi- 
ciency in destroying a np condensate. In such a situation, 
despite the fact that the isospin asymmetry may be very 
large, a np condensate survives and exists in the form of 
a Bose-Einstein condensate of deuterons. 

The transition from BCS superconductivity to Bose- 
Einstein condensation (BEC) occurs in a Fermi system, 
if either density is decreased or the attractive interaction 
between fermions is increased sufficiently. This transi- 
tion was studied, first, in superconducting semiconduc- 
tors Q and then in an attractive Fermi gas |9|. Later 
it was realized that an analogous phase transition takes 
place in symmetric nuclear matter, when np Cooper pairs 
at higher densities go over to Bose-Einstein condensate 
of deuterons at lower densities 0, |3 ■ During this tran- 
sition the chemical potential changes its sign at certain 
critical density (Mott transition) , approaching half of the 
deuteron binding energy at ultra low densities. In Ref. |5j 
crossover from np superfluidity to BEC of deuterons was 
investigated in the T-matrix approach, where the pole 
in the T-matrix determines the critical temperature of 
BEC of bound states in the case of negative chemical po- 
tential and the critical temperature of the appearance of 
np Cooper pairs in the case of positive chemical poten- 
tial. The influence of isospin asymmetry on the BEC- 



BCS crossover in nuclear matter was studied in Ref. [l(]j 
within the BCS formalism. It has been shown that Bose- 
Einstein condensate is weakly affected by an additional 
gas of free neutrons even at very large asymmetries. The 
same conclusion was also confirmed in Ref. |llj on the 
base of the variational approach for the thermodynamic 
potential. 

The recent upsurge of interest to the BEC-BCS 
crossover is caused by finding the B CS p airing in ul- 
tracold trapped quantum atom gases \\A |lj|. In this 
study we examine the possible signatures of the BEC- 
BCS crossover in low-density nuclear matter. It may have 
interesting consequences, for example, in the far tails of 
the density profiles of exotic nuclei, where a deuteron con- 
densate can exist in spite of the fact that the density there 
can be quite asymmetric. Besides, similar physical effects 
can play an important role in expanding nuclear matter, 
formed in heavy ion collisions, or in nuclear matter in 
the crust of a neutron star. The main emphasis is laid on 
the behavior of the density, spin and isospin correlation 
functions across the BEC-BCS transition region. The 
study is motivated by the results of Ref. fl^l, where the 
authors state that the density correlation function of a 
two-component ultracold fermionic gas of atoms changes 
sign at low momentum transfer and this represents an 
unambiguous signature of the BEC-BCS crossover. This 
statement is checked for nuclear matter taking into ac- 
count additional factors: finite isospin asymmetry or fi- 
nite temperature. In both cases, this criterion fails to 
provide a correct description of the density-driven BEC- 
BCS crossover and cannot serve as the universal feature 
of transition between two states of the system. 

Basic equations. Superfluid states of nuclear matter 
are described by the normal / and anomalous g distribu- 
tion functions of nuclcons 

/kiK2 = Tr£ia+ Cki, .g^iKa = TrgiaK^a^i, (1) 

where k = (k, ct, r), k is momentum, <t{t) is the projec- 



tion of spin (isospin) on the third axis, g is the density 
matrix of the system. We shall study np pairing cor- 
relations in the pairing channel with total spin S and 
isospin T of a pair S ~ 1, T ~ Q and the projections 
Sz = Tz = Q. In this case the distribution functions for 
isospin asymmetric nuclear matter have the structure 



/(k) = /oo(k)croTo + /o3(k)croT-3, 

5(k) =53o(k)cr30-2r2, 



(2) 



where Ui and Tk are the Pauli matrices in spin and isospin 
spaces, respectively. Using the minimum principle of the 
thermodynamic potential and procedure of block diago- 
nalization [3, one can obtain expressions for the distri- 
bution functions 



/oo(k)4-J-(tanh§ + tanhf), 

fosik) = - ( tanh -^ - tanh -^ ) , 

A(k) / E+ E- \ 

53o(k) = -^(tanh^+tanh^). 



Here 
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Ek ± Sfi - ^^[+AHk)±Sfi, ^k 
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(3) 
(4) 
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M,(6) 



A being the energy gap in the quasiparticle excitation 
spectrum, m being the effective nucleon mass, /i and dfi 
being half of a sum and half of a difference of neutron 
and proton chemical potentials, respectively. 

Equations, governing np pairing correlations in the 
S — 1, T = pairing channel, can be obtained on 
the base of Green's function formalism and have the 
form dH [13 

k' '^' 

(7) 



Ek' 



^9=yY.[f^Ek)-f{Et) 



(8) 
(9) 



where f{E) is Fermi distribution function. Eq. Q is 
equation for the energy gap A and Eqs. (jHl), © are equa- 
tions for the total density g ~ gp + gn and neutron excess 
6g = gn — gp = ag {a being the asymmetry param- 
eter). Note that, since we consider unitary superfluid 
state (AA+ ex /), Eqs. 10-® formally coincide with the 
equations for a two-component isospin asymmetric super- 
fluid with singlet spin pairing between unlike fermions. 
Introducing the anomalous density 



^(k) =< 
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A(k)/ 
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and using Eq. ((SJ, one can represent Eq. Q for the en- 
ergy gap in the form 



m 



ij{k) + {1- uk) J2 ^(k, kOV^k') = 2A.7^(k). (10) 



In the limit of vanishing density, Uk -^ 0, Eq. H10() 
goes over into the Schrodinger equation for the deuteron 
bound state |j,|lOj. The corresponding energy eigenvalue 
is equal to 2iJ,. The change in the sign of the mean chemi- 
cal potential fi of neutrons and protons under decreasing 
density of nuclear matter signals the transition from the 
regime of large overlapping np Cooper pairs to the regime 
of non-overlapping bound states (deuterons). 

Let us consider the two-body density correlation func- 
tion 

P(x,x') = Tr£iA?T,(x)Ar7.(x'), An(x) = n(x) - h, (11) 

err 

1 v^ 
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Its general structure in the spatially uniform and 
isotropic case reads 16] 



P(x, x') = g8{v) + gD(r), r = x - x' 



(12) 



The function D(r) is called the density correlation func- 
tion as well. We will be just interested in the behavior of 
the function -D(r). The trace in Eq. Hl()|l can be calcu- 
lated, using definitions (Q and Wick rules. Taking into 
account Eqs. |5J and going to the Fourier representation 



0{q) 



d^re'^^^Dir), 



one can get 



D{q)=lf>{q)-lf>{q)-l'/{q), 



(13) 



where 



^/°(9) ^^J,^drr^jo{rq)[j^dkefoo{k)Mrk) 
^f (9) ^^jrdrr^Joirq)[!,^dkk'fo3ik)Mrk) 
^f (9) = ^ /o°° dr r^Mrq) [/^ dk eg,,{k)j,{rk) 



Here jo is the spherical Bcssel function of the first kind 
and zeroth order. The functions /9'',/9'^ and l¥^ repre- 
sent the normal and anomalous contributions to the den- 
sity correlation function. Analogously, we can consider 



the two-body spin correlation function 
5^i.(x,x') = TrgiAsp(x)As,,(x'), As^(x) = Sp(x) - s^, 

(14) 

crcr'r 

- — V e^('''"*')''a+ (a ) ^aw / 
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and the two-body isospin correlation function 
r^,(x,x') = TreAt^(x)Ai(x'), AtV(x) = i^{x) - t^, 

(15) 

*"mW = 2 XI ^<^r(x)(T,Orr'V'<Tr'(x) 
crrr' 

~ 2V^ ^ "kcTTVMJT-T'«k'<TT', 

frrr'kk' 
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Their general structure for isospin asymmetric nuclear 
matter without spin polarization is 

S^A^,^')^^S^,S{r) + gS^Ar), (16) 

r^,(x,x') = I S^, S{r) + ^ ie,,,3 -^(r) + eT,,,(r). (17) 

Then, calculating traces in Eqs. p4|l . p5|l . for the Fourier 
transforms of the spin and isospin correlation functions 
one can get 

5^.(9) = -i{V(^/"(?)+^f('Z)) (18) 






(19) 



Note that if to put i^ = ^ = 3 in Eqs. (jIHl), lO, one 
gets the longitudinal spin 5' and isospin T' correlation 
functions, while setting fi,!/ — 1,2 gives the transverse 
spin and isospin correlation functions 

SIM = -%(^r(9) + iTii) + ^f (9)) ^ V5*(g), 



//,!/= 1,2, 

(20) 
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KAl) = -f{lT(l)~ir(l)+l',"i'i)) - S^'^T^q). 



The following relationships between the correlation func- 
tions hold true 

S\q)^^, S\q)^T\q). (21) 

At zero temperature and zero momentum transfer, the 
correlation functions satisfy the sum rule 



S\q = Q)=T'{q = 0) 
1 



(22) 



27r2g 



dkk\f^oik) + f^3ik)+gio{k)) 



where the r.h.s. is independent of density and isospin 
asymmetry. Besides, the transverse isospin correlation 
function satisfies the relationship 



T\q = 0) = - 



1 



2tt^q 

1-Q 



dkk^fSoik) - fS,{k) + gloik)) 

(23) 



where the r.h.s. is independent of density. 

Correlation functions in nuclear matter with 
a np condensate. Further for numerical calculations 
we shall use the effective zero range force, developed in 
Ref. [131 to reproduce the pairing gap in 5 = 1,T = 
pairing channel with Paris NN potential: 



y(ri,r2) =uo 



{!-,( 



Qi 



ri+r2 



go 



> 



I"! -I"2), 



(24) 



where qq ^ 0.16 fm^ is the nuclear saturation density, 
Vq = — 530MeV • fm , 77 = 0, tti = mg, niQ being the 
effective mass, corresponding to the Gogny force DIS. 
Besides, in the gap equation Q, Eq. (|^ must be sup- 
plemented with a cut-off parameter, Ec = 60 MeV. Just 
this set of parameters, among total three parametriza- 
tions, used in Ref. [l3|, corresponds to the formation of 
bound states at nonzero energy in low-density region of 
nuclear matter. 

To find the correlation functions one should first solve 
the gap equation ([TJ self-consistently with Eqs. (jSj), lO. 
Then the correlation functions can be determined di- 
rectly from Eqs. ^^, lO and ^^. The results of nu- 
merical determination of the energy gap as a function of 
density for different asymmetries at zero temperature are 
shown in Fig. 1. As one can see, with increasing asymme- 
try the magnitude of the energy gap is decreased and the 
density interval, where a np condensate exists, shrinks to 
lower density. In reality solutions exist for any a < 1 
(the phase curves for larger values of a are not shown 
in Fig. 1) and correspond to the formation of BEG of 
deuterons at very low densities of nuclear matter. 

Now we consider the correlation functions D{q) and 
S^{q) for symmetric nuclear matter at zero temperature, 
depicted in Fig. 2 (at a = 0, T\q) = S\q)). The den- 
sity correlation function changes sign at low momentum 
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FIG. 1: Energy gap as a function of density at zero tempera- 
ture and different asymmetries. 



transfer when the system smoothly evolves from the BEC 
regime to the BCS one. These two regimes are distin- 
guished by negative and positive values of the chemical 
potential [i, respectively. In view of Eq. if^ , the longitu- 
dinal spin correlation function S"' (q) changes sign through 
the BEC-BCS crossover as well. The transverse spin cor- 
relation function, and, according to Eq. (|21|) . the longitu- 
dinal and transverse isospin correlation functions change 
fluently between BEC and BCS limits. The behavior 
of the density correlation function in isospin symmetric 
case at zero temperature qualitatively agrees with the be- 
havior of the density correlation function in an ultracold 
fermionic atom gas with singlet pairing of fermions [l5j . 
In Ref. [l5j , the change in the sign of the density correla- 
tion function at low momentum transfer was considered 
as a signature of the BEC-BCS crossover. We would like 
to extend their calculations taking into account the finite 
isospin asymmetry and finite temperature. 

Fig. 3 shows the dependence of the density correlation 
function Z?((7 = 0) at zero momentum transfer as a func- 
tion of density for a set of various isospin asymmetry 
parameters and zero temperature. It is seen that with 
increasing the asymmetry parameter the density corre- 
lation function decreases. For strong enough asymme- 
try, the function D{ci = 0) is always negative. In ac- 
cordance with the above criterion, the density region, 
where the function £'(g = 0) has positive or negative 
values, would correspond to the BEC or BCS regime, 
respectively. Hence, as follows from Fig. 3, for strong 
isospin asymmetry we would have only the BCS state 
for all densities where a np condensate exists. Obviously, 
this conclusion contradicts with the behavior of the mean 
chemical potential /i, being negative at very low densi- 
ties for any a < 1, and, hence, giving evidence to the 
formation of BEC of bound states [llj. Thus, at strong 
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FIG. 2: Density and transverse spin correlation functions as 
functions of momentum at zero temperature and different 
densities for symmetric nuclear matter. 




0.00 0.04 0.08 0.12 0.16 

D rfm'l 



0.20 



FIG. 3: Density correlation function D(q = 0) as a function 
of density at zero temperature for different isospin asymmetry 
parameters. 



enough isospin asymmetry the criterion of the crossover, 
based on the change of the sign of the density correla- 
tion function, fails to predict the transition to the BEC 
of deuterons in low-density nuclear matter. 

Now we consider symmetric nuclear matter at finite 
temperature. Fig. 4 shows the dependence of the den- 
sity correlation function D{q = 0) at zero momentum 
transfer as a function of density for a set of various tem- 
peratures. It is seen that for not too high temperatures 
the density response function is nonmonotonic and twice 
changes sign in the region of low densities. Hence, in 
accordance with the above criterion, we would have the 
density interval Qi < Q < Q2 with the BEC state, sur- 
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FIG. 4: Density correlation function D(q = 0) as a function of 
density at different temperatures for symmetric nuclear mat- 
ter. 



ties of fermions of different species) or finite tempera- 
ture. Hence, the criterion in Ref. 15] cannot be consid- 
ered as the universal indication of the BEC-BCS tran- 
sition. During the Mott transition, when the chemical 
potential changes sign, there is a qualitative change in 
the quasiparticle energy spectrum: the minimum shifts 
from a finite (BCS state) to zero-momentum value (BEC 
state) (see Eq. (0 and Ref. 1,7]). As such, the presence 
(BCS) or absence (BEC) of the singularity in the mo- 
mentum distribution of the quasiparticle density of states 
represents the universal signature of the BEC-BCS tran- 
sition. This transition may be relevant, and could give a 
valuable information on jip pairing correlations, in low- 
density nuclear systems, such as tails of nuclear density 
distributions in exotic nuclei, produced at radioactive nu- 
clear beam facilities, expanding nuclear matter in heavy 
ion collisions, low-density nuclear matter in outer regions 
of neutron stars, etc. 



rounded by the density regions with the BCS state. How- 
ever, this conclusion contradicts with the behavior of the 
mean chemical potential fj, for these temperatures, being 
a monotone function of density and indicating the for- 
mation of a BEC state at low densities (^ < 0) and a 
BCS state at larger densities (/i > 0). Thus, at finite 
temperature the criterion of the crossover, formulated in 
Ref. 15], fails to provide the correct description of the 
transition between two regimes. 

In summary, we have calculated the density, spin and 
isospin correlation functions in superfluid nuclear matter 
with np pairing correlations, intending to find the pos- 
sible signatures of the BEC-BCS crossover. It is shown 
that the transverse spin, and longitudinal and transverse 
isospin correlation functions satisfy the sum rule at zero 
momentum transfer and zero temperature, and change 
smoothly between BEC and BCS regimes. In Ref. [l5j, 
it was learned that the density correlation function in a 
two-component ultracold fermionic atom gas with singlet 
pairing of fermions changes sign at low momentum trans- 
fer across the BEC-BCS transition, driven by a change in 
the scattering length of the interaction at zero temper- 
ature. We have shown that for spin triplet pairing the 
longitudinal spin correlation function plays an analogous 
role to the density correlation function and changes sign 
at low momentum transfer across the crossover in sym- 
metric nuclear matter at zero temperature. However, 
while giving a satisfactory description of the density- 
driven BEC-BCS crossover in dilute nuclear matter at 
zero temperature for the isospin symmetric case, this 
criterion fails to provide the correct description of the 
crossover at finite isospin asymmetry (nonequal densi- 
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